Momentum-resolved single-particle spectral function for TiOCl from a combination of 
density functional and variational cluster calculations 



o 
o 

< 

00 
(N 



c 

o 

CM 
> 
m 
m 

d 

00 

o 



M. Aichhorn, 1 ' 2 T. Saha-Dasgupta, 3 R. Valentf, 4 S. Glawion, 5 M. Sing, 5 and R. Claessen 5 

J Insitut fur Theoretische Physik und Astrophysik, 
Universitat Wiirzburp, Am Hubland, 97074 Wiirzburg, Germany 
2 Centre de Physique Theorique, Ecole Poly technique, CNRS, 91128 Palaiseau Cedex, France 
3 S.N. Bose National Centre for Basic Sciences, JD Block, 
Sector III, Salt Lake City, Kolkata 700098, India 
Jf Institut fur Theoretische Physik, Goethe- Universitat Frankfurt, 60438 Frankfurt/Main, Germany 
5 Experimentelle Physik 4, Universitat Wiirzburg, Am Hubland, 97074 Wiirzburg, Germany 

We present results for the momentum-resolved single-particle spectral function of the low- 
dimensional system TiOCl in the insulating state, obtained by a combination of ab initio Density 
Functional Theory (DFT) and Variational Cluster (VCA) calculations. This approach allows to 
combine a realistic band structure and a thorough treatment of the strong correlations. We show 
that it is important to include a realistic two-dimensional band structure of TiOCl into the effective 
strongly-correlated models in order to explain the spectral weight behavior seen in angle-resolved 
photoemission (ARPES) experiments. In particular, we observe that the effect of the interchain 
couplings is a considerable redistribution of the spectral weight around the V point from higher 
to lower binding energies as compared to a purely one-dimensional model treatment. Hence, our 
results support a description of TiOCl as a two-dimensional compound with strong anisotropy and 
also set a benchmark on the spectral features of correlated coupled-chain systems. 

PACS numbers: 71.27.+a,71.10.-w,71.10.Fd 



I. INTRODUCTION 

In recent years a significant amount of research has 
been dedicated to strongly-correlated materials with re- 
duced dimensionality since they exhibit a large variety 
of fascinating dimension-related properties. An exam- 
ple is the layered quantum spin system TiOCl, where 
bilayers of Ti-0 are separated by Cl~ ions. This sys- 
tem was originally thought to be a possible candidate 
for a RVB superconductor upon doping^ because of its 
frustrated triangular lattice geometry. Later on, various 
experimental measurements^ 3 -'^' 5 -^ revealed that TiOCl 
shows in fact an anomalous spin-Peierls behavior with 
two consecutive phase transitions. Magnetic susceptibil- 
ity was initially described in terms of a one-dimensional 
spin-1/2 Heisenberg model with a large intra-chain cou- 
pling constant J w 700 K^. It is though well known 
that susceptibility is not very sensitive to different mod- 
els and recent ab initio DFT studies^ showed that the un- 
derlying interactions for this system can be understood 
in terms of a spin-1/2 Heisenberg model with strong in- 
trachain antifcrromagnctic interactions J\ — 6607^ and 
weaker interchain ferromagnetic interactions J2 = — 16K, 
J3 = — 10K. This model reproduces the magnetic sus- 
ceptibility measurements and sets a framework for un- 
derstanding the puzzling spin-Peierls phase transitions in 
this compound. Only recently, research has also focused 
on high-pressure studie s 10 ' 11 ' 12 ' 13 as a possible way to 
drive the system metallic. 

At room temperature and ambient pressure, the 
system is a Mott insulator with a charge gap of 
about 2eV&^- The electronic structure in this high- 
temperature phase has been examined by angle-resolved 



photoemission spectroscopy (ARPES) . 14 ' 15 In agreement 
with previous experimental evidence, the results show a 
strong anisotropy of the correlated band structure, with 
significant dispersion of the Ti 3d bands along the chains 
(crystallographic 6-direction) , and almost flat bands per- 
pendicular to the chains. 

On the theoretical side, the electronic properties of Ti- 
OCl have been studied by means of ab initio DFT cal- 
culations within the local density approximation (LDA) , 
the LDA+U^, B3LYR±£ and also in combination with 



the dynamical mean-field theory (DMFT) 
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which 



is a modern method for dealing with strong correlations. 
It was shown that a proper treatment of non-local corre- 
lations is crucial for a reasonable description of the single- 
particle gap^i 

However, the momentum dependence of the spectral 
function A(k, u>) seen in ARPES is still puzzling. It has 
been shown that an ab initio calculation without proper 
treatment of correlations is insufficient i^ ' 15 ' 16 ' 17 On the 
other hand, describing the compound by a simplified 
one-dimensional strongly correlated model was not suc- 
cessful either*^ Furthermore, LDA+DMFT could so far 
only produce the momentum integrated local density of 
states (DOS) without any information on the momentum 
dependence of the spcctrai 18 ' 19 ' 20 This situation, having 
no calculation for the momentum-resolved spectral func- 
tion A(k, to) at hand, is partly due to the fact that there 
are only few methods that can deal with all the require- 
ments of such calculations. This work is intended to fill 
this gap and investigates the influence of the true two- 
dimensional band-structure on the momentum-resolved 
A(k, oj) in the presence of strong correlations. A success- 
ful technique for this purpose is the Variational Cluster 
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Approach (VCA)2L22. 

In what follows we apply a two-step procedure to study 
the spectral function, as has been proposed by Chion- 
ccl et First. DFT calculations within the LDA are 

carried out, and localized Wannier functions are con- 
structed by the N-th order muffin-tin-orbital (NMTO)^ 
downfolding technique. Using the LDA Hamiltonian ex- 
pressed in these Wannier functions as the non-interacting 
part, and adding Coulomb and Hund interaction terms, 
we arrive at the correlated low-energy model. By ap- 
plying VCA to this model Hamiltonian, we show that 
the inclusion of the inter-chain processes leads to a sig- 
nificant redistribution of spectral weight from higher to 
lower binding energies. Since these processes enhance 
the asymmetry of the strongly-correlated band structure, 
they are crucial for the reproduction of the asymmetric 
bands seen in ARPES measurements. Our calculations 
show that the Hubbard-modcl description is appropriate 
for TiOCl if effects beyond the one-dimensional descrip- 
tion are included. Moreover these results should be valid 
for a large variety of correlated low-dimensional coupled- 
chain systems. 

The paper is organized as follows: In Scct.|TT]wc discuss 
the construction of the low-energy Hamiltonian, as well 
as the VCA, which is subsequently used for the calcula- 
tion of the correlated spectral function. Sect. lIIII contains 
our results of the multi-band as well as of the single-band 
Hubbard model and in Sect. IIVI we present our discus- 
sions and conclusions. 



II. THEORY 

In many transition metal oxides electronic correlation 
effects are very important for a proper description of the 
physical properties. However, it is a known fact that first- 
principle calculations suffer from an insufficient treat- 
ment of these effects. In order to take the strong cor- 
relations into account in our calculation for TiOCl, we 
apply a two-step procedure (LDA + VCA) that has first 
been introduced by L. Chionccl et al.£2- It consists of 
the construction of the correlated low-energy Hamilto- 
nian based on density-functional theory on the one hand, 
and the solution of the resulting low-energy Hamilto- 
nian using the VCA on the other hand. In Ref. H^, the 
authors study the non-quasiparticle states in the half- 
metallic compound Cr02 and find good agreement with 
experiments. Moreover, a comparison with LDA+DMFT 
calculations showed the applicability of the LDA+VCA 
approach. Recently, it has also been used to explain the 
pscudogap in TiN, where also the momentum-resolved 
spectral function has been calculated^ 



A. Low-energy Hamiltonian 

For a complete description of the electronic structure 
of a given material it is in principle necessary to consider 



all electronic degrees of freedom of the underlying con- 
stituents. Calculations within DFT can to some extent 
fulfill this requirement. However, it is clear that only 
certain states and orbitals contribute to the low-energy 
physics. For this reason one is interested in finding an 
effective model that describes the low-energy physics on 
the one hand sufficiently accurate and has, on the other 
hand, not too many degrees of freedom. 

In the present case of TiOCl, DFT calculations within 
the LDA approximation have shown that the relevant or- 
bitals at low energies are the Ti 3d orbitals, which are 
split into tig and e g manifolds due to the octahedral crys- 
tal field provided by the ligands. Since the Ti 3+ ion is 
in a 3d 1 configuration, the relevant states closest to the 
Fermi energy are of predominantly tig character. 

For the construction of the low-energy Hamiltonian, 
we performed DFT calculations within the LDA using 
the linearized muffin-tin-orbitals (LMTO) basis set. The 
localized orbitals, which are the basis of the interacting 
model, are constructed using the NMTO method. By 
using the downfolding technique^, the hybridization of 
the Ti-tig orbitals with the ligand orbitals (O-p and Cl-p) 
are taken into account, yielding an effective set of ti g or- 
bitals. These orbitals represent the LDA band structure 
with great accuracy, and are used as the non-interacting 
part of the many-body Hamiltonian. The matrix ele- 
ments of the NMTO Hamiltonian # LDA (k) in the basis 
set of localized NMTO Wannier functions give the trans- 
fer integrals , and the non-interacting Hamiltonian 
can be written as 

The indices label the lattice sites by i, j, as well as the 
tig orbitals by a, (5, and tr denotes the spin. 

To include correlation effects into the low-energy de- 
scription, we add interaction terms to the Hamiltonian, 

H =Hq — [l ^ ' Ilia H — ^ ] n iaa n iaa 
ia icta 

+ Y n ia n i(3- JZ S L S h (2) 

For convenience, we introduced the chemical potential fi 
in the Hamiltonian. We will refer to this Hamiltonian 
as tig model, and give all energies throughout the pa- 
per in units of electron volt (eV). The full low-energy 
model Eq. @ consists of the single-particle terms (Ho 
and /x), the diagonal (density-density) interactions (U, 
U', and J z ), and the non-diagonal (spin- flip) term (J, 
third line). In this study we consider only the case 
J z = J. rii a = riia^ + riiai is the orbital occupation 
operator, and S? a , S^~ a , S~ a are the components of the 
spin-i operator on site i in orbital a. The interaction 
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parameters U, U' , and J are not independent, but fulfill 
the relation U' = U — 2J. At this point it is impor- 
tant to note that we can include a full SU (2) symmetric 
exchange term. Since the method we consider is not af- 
fected by any sign problem, it has no restriction on the 
type of couplings that can be included.— For the inter- 
action parameters U and J one can find several values 
in the literature, ranging from U = 3.0 cV to U — 4.0 eV 
and J = 0.5 cV to J = l.OeV^MO&l&SO Since we want 
to study also the influence of these parameters on the 
single-particle properties, we have performed calculations 
with different values, and indicate the actual value at the 
corresponding location in the paper. 

In this work, we also address the question whether the 
orbital degrees of freedom are important for the low- 
energy physics or not. Since TiOCl does not crystal- 
lize in a perfect cubic symmetry, the threefold degener- 
acy of the tig manifold is lifted. LDA+U 2 -^ and also 
LDA+DMFT— calculations have shown that the ground 
state shows predominantly d xy character (the local refer- 
ence frame is z = a, and x and y axes rotated by 45° with 
respect to b and c) , with only very small admixture of the 
other orbital degrees of freedom, a picture that we will 
also find in our following calculations. This is in contrast 
to IPT-DMFT calculations where a sizable admixture 
of the other orbital degrees of freedom is found. 

In order to investigate the effective one-band model 
that consists of the d xy orbital only, we performed a 
NMTO downfolding procedure integrating out all other 
degrees of freedom, and keeping only the d xy channel. 
In this one-band model, the only interaction terms are 
the ones proportional to the Hubbard onsite U, and the 
low-energy one-band Hamiltonian finally reads 



(3) 



where the Hubbard interaction U is the same as for the 
tig -model. The effective hopping parameters Uj arc 
again the matrix elements of £T LDA (k) in the Wannier 
basis set. 
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FIG. 1: Triangular lattice structure and two possible clusters 
tiling the lattice, i) 4-site cluster including inter-chain self- 
energies, ii) ID clusters neglecting inter-chain self-energies. 
Full circles and solid lines mark sites and bonds inside a clus- 
ters. Next-next-nearest neighbor hopping tz is only drawn 
once for clarity. In the case of the ta g manifold each lattice 
site consists of three orbitals. 



rection. It gives just a constant shift in energy which can 
be absorbed in the chemical potential. 

As mentioned in the introduction, the VCA is a quan- 
tum cluster method capable of treating strong short- 
ranged correlations. The main idea is to approximate the 
sclf-cncrgy of the original model, which is defined on an 
infinite lattice, by the self-energy of a finite cluster, the 
reference system. The variational principle states that 
the optimal solution is given by the stationary points of 
the grand potential fi[E] as a function of the self-energy 
S. Parametrizing the self-energy by the single-particle 
parameters t' of the reference system, we can write the 
grand potential as 

fi(t') = fi' + Trln(G^ - S(t')) -1 - Tr In Gf , (4) 

where f2' and Gf are respectively the grand canonical 
potential and the Green's function of the reference system 
and Go,t is the non-interacting Green's function of the 
physical (lattice) system. The stationary condition reads 



B. Variational Cluster Approach 
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After having constructed the low-energy Hamiltonian 
using ab initio techniques, we use the VC A 21 ' 22 in order 
to calculate the spectral function of this model. Since 
we deal with an effective low-energy Hamiltonian that 
involves no other uncorrclatcd ligand states (O-p, Cl-p), 
but only the correlated Ti-tig orbitals, the application of 
VCA is straightforward, and is from a technical point of 
view exactly equivalent to standard multi-orbital calcu- 
lations for Hubbard- model Hamiltonians. The only dif- 
ference is that the non-interacting part is determined by 
the procedure discussed in the previous subsection. Fur- 
thermore, since there are no explicit ligand states in the 
Hamiltonian, there is no need for a double-counting cor- 



It is important to note that the interaction parameters 
are not variational parameters, since, by construction of 
the VCA, the interaction terms of the reference system 
and the original lattice model must not differ. In this 
study, we restrict ourselves to local interactions only since 
the VCA in its strict sense cannot be used for models with 
non-local interactions without further approximations. 

In its spirit, the VCA is closely related to the dynam- 
ical mean-field theory (DMFT), where in the latter case 
the self-energy is obtained from an impurity problem. 

The actual VCA calculation is done in the following 
steps. First, we determine the ground state of the ref- 
erence system, i.e. a cluster of finite size as depicted in 



4 



Fig. [T] The interacting Green's function is calculated, 
and since the non-interacting Green's function of the ref- 
erence system is known, the self-energy can readily be 
obtained using Dyson's equation. Using the grand po- 
tential f2' of the reference system, the Green's function 
Gf, and the self-energy S(t')), Eq. [4] is evaluated us- 
ing the technique of Q-matrices<2I Note that the Green's 
functions Go,t, Gf and the self-energy S(t') in Eq. ^ 
are matrices not only in site and spin indices, but also 
carry an orbital index. In fact, this is the only difference 
of the application of VCA in the present case compared 
to the numerous previous applications to the single-band 
Hubbard model. 

As a reference system solver, we use the Band-Lanczos 
exact diagonalization technique at zero temperature, 
which means that for the full ti g manifold, we can eas- 
ily consider clusters with at most 4 sites, yielding a 12- 
orbital Green's function Gf . For the single-band model, 
we consider clusters up to 12 sites. We exploit parti- 
cle number and spin conservation, therefore the sizes of 
the largest Hilbert spaces that we have to consider are 
N = 14520 states in the 4-site multi-orbital case, and 
TV = 853776 states in the 12-sitc single-orbital case, re- 
spectively. Since we are considering an exact diagonal- 
ization method for solving the cluster problem, all in- 
teractions in the Hamiltonians Eq. ((5J) and Eq. ([3]) are 
treated exactly and on the same footing. This is a clear 
advantage compared to, e.g., using the Hirsch-Fye quan- 
tum Monte-Carlo method as impurity solver, since in the 
latter case approximations to the interaction terms of the 
Hamiltonian have to be done.— 

The VCA approach has been tested thoroughly and 
used successfully for many investigations in recent years. 
Several studies on the cuprate-based high-temperature 
superconductors have shown that this approach can re- 
produce salient features of these materials, such as the 
ground-state phase diagram , 27 i 28 ' 29 i 30 or the opening of 
the pseudogap at low hole doping, accompanied with the 
occurrence of Fermi arcs , 30 ' 31 in very good agreement 
with experiments and results obtained by the cellular dy- 
namical mean-field theory (CDMFT) (see, e.g, Refs. l32l). 
Recently, the VCA could also reproduce the pairing sym- 
metry of the iron-based superconductors.— 

The VCA has also been used for multi-orbital sys- 
tems, which is relevant for the combination with ab- 
initio methods. On the pure methodological level, 
the metal-insulator transition in infinite dimensions was 
studie d 34 ' 35 , and very good agreement with dynamical 
mean-field calculations was found. An application to 
real materials was done in Refs. l36l|37l where the com- 
pounds NiO, CoO, and MnO have been studied and very 
good agreement with experimental photo-emission data 
has been found. 

Details on the practical implementation of the VCA, 
inclu ding t ests and benchmarking, can be found in 
Refs. I2i|27|3il3ili31il 

In general, all the single-particle parameters t' are vari- 
ational parameters of the VCA. In practice, one chooses 



a physically motivated subset in order to keep the nu- 
merical calculations feasible. Here we make the following 
choice. For a thcrmodynamically consistent description 
of the densities, it is crucial to consider the onsite ener- 
gies, i.e. the local terms of the single-particle Hamilto- 
nian e' a = (tf")', as variational parameters.— We define 
the average e' = \{s' xy + £ yz ) and the crystal-field split- 
ting A' cf = e' yz — e' xy , which are then used as the varia- 
tional parameters of the VCA. Note that in the single- 
band case, Eq. ([3]), one has to deal with e' only. One 
has to be aware that the variational parameter A' ci does 
not impose an artificial orbital polarization of the system, 
since it is a parameter in the variational procedure and no 
physical external field. Hence, using e' and A' c{ , the or- 
bital occupancies are determined in a fully self-consistent 
way. 

The main property investigated in this work is the 
single-particle spectral function which we define as 

A(k,uj) = — trlmG(k,u;) , (6) 

Since we broke the translational invariance of the sys- 
tem by introducing the cluster tiling, a proper pcriodiza- 
tion of the lattice quantities is needed in order to restore 
translational symmetry, an issue also important in cluster 
DMFT calculations 42- Here, we choose to use the peri- 
odization of the Green's function since the periodization 
of the self-energy gives unphysical results in the insulat- 
ing phased In other words, starting from the Green's 
function that depends on two momenta, G(k, k',w), one 
restores the fully translationally invariant Green's func- 
tion G(k, ui), by neglecting the off-diagonal elements, and 
taking k = k' only. It has been shown that this is a well 
justified approximation to calculating the momentum- 
dependent spectral function^ 3 - The Green's function G 
is in general a matrix in orbital indices and A(k, lu) is 
given by the trace over the orbital degrees of freedom. 



III. RESULTS 
A. Full t-zg model vs. effective one-band model 

Before we come to a detailed analysis of the single- 
band model, we first want to check the validity of the 
restriction to the lowest d-orbital. 

For this reason, we performed ab initio calculations to 
determine the full single-particle Hamiltonian of the sys- 
tem. Since Ti 3+ is in a 3d 1 configuration, the e g orbitals 
are unoccupied and can be projected out, and the full ki- 
netic Hamiltonian is downfolded to the threefold degener- 
ate t2 g manifold. Results show, that the system exhibits 
a strong anisotropy, with the largest hopping integrals in 
the crystallographic 6-direction, see Fig.[IJ almost one or- 
der of magnitude larger than the other transfer integrals. 
Moreover, the threefold degeneracy is lifted and the man- 
ifold split into the lower d xy and the higher d xz and d yz 
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FIG. 2: Comparison of the single-particle spectral function 
of the t2 g model (top) and the single-band model (bottom), 
both calculated with a 2 x 2 reference system, see Fig. [T] i) . 
Parameters are U = 3.3 eV, J = 0.5 eV. For the hopping 
parameters see text. The chemical potential has been chosen 
such that i) the system is insulating with n = 1.0, and ii) the 
position of the occupied states coincide in both calculations. 
Lorentzian broadening of r\ = 0.02 eV was used. 

orbitals. Note that for the orbital designation we con- 
sider the same local reference frame as in Ref. Il4lfl8l . The 
crystal-field splitting between the ground state and the 
first excited state obtained from LDA is about 0.42 eV. 
This theoretical value is in reasonable agreement with ex- 
perimental results^i^ Despite this splitting, the orbital 
sector is not fully polarized in the LDA calculations, and 
the occupation of d xy is about 0.49, with 0.51 electrons 
in the other two orbitals. 

In order to perform our LDA+VCA calculations, we 
take the downfolded Hamiltonian of the ab initio calcu- 
lations, and add the interaction and exchange terms ac- 
cording to Eq. In the upper panel of Fig. [2] we show 
the results for the spectral function, Eq. ([6]), calculated 
for this three-band t2 g model using typical parameters 
U = 3.3 eV and J = 0.5 eV. The bands which are located 
just above the Fermi level, between roughly 0.5 eV and 
2.0 eV, have d xz and d yz character, and remain almost 
unchanged by the strong interactions. 

The behavior of the d xy orbital is strikingly different. 
It splits into two bands that can be identified with the 
lower and upper Hubbard band, located roughly around 
— 1.0 eV and 2.5 eV, respectively. By inspecting the terms 
of the Hamiltonian related to the crystal-field splitting, 
i.e. H c { = A c f J2i( n yz + n zx — n xy ), we can calculate the 
orbital polarization p = dfl/dA c f. We find a value of 
p = —0.99, meaning that the system is almost perfectly 
polarized into the d xy orbital, which is in agreement with 
recent LDA+CDMFT calculations^ Interestingly, this 
polarization is found without any sizable increase of the 



FIG. 3: (Color online) Comparison of the spectral function of 
the t2 g model for two different sets of interaction parameters. 
Top: U = 4.0 eV, J = 0.7 eV. Bottom: U = 4.0 eV, J = 
0.5 eV. The vertical dashed lines mark the edges of the single- 
particle gap. Lorentzian broadening of n = 0.02 eV was used. 



crystal- field splitting in the variational procedure, i.e., 
A c f w A' ci , but it is only due to the inclusion of strong 
local interactions. 

This result gives rise to the question, to which extent 
a single-band Hubbard model can describe the occupied 
states relevant for comparison with ARPES. We took pa- 
rameters from a full downfolding to the Ti d xy orbital 
only, cf. first column of Table I in Ref. [H. Using the 
same value of U = 3.3 eV we calculate the spectral func- 
tion for Hamiltonian Eq. ([3]). The results are shown in 
the lower panel of Fig. [5J In order to avoid effects coming 
from different cluster sizes, we used also a 2 x 2 cluster 
for this comparison. Note that below the Fermi level the 
agreement between the single-band model and the d xy 
part of the t2 g model is excellent. For this reason we 
conclude that for a comparison of spectra with experi- 
mental ARPES measurements the Hamiltonian Eq. Q 
is a reasonable starting point. 

Before turning to a more detailed analysis of the spec- 
tra obtained from Eq. ([3]), let us briefly discuss the single- 
particle gap A, defined as the energy difference between 
the lowest unoccupied and the highest occupied state. 
For a comparison of this quantity with experiments, it is 
clear that the single-band model is not sufficient, since it 
does not describe the excited states in the ti g manifold. 
However, we extracted the gap A from the spectral func- 
tion of the t2 g model, for different values of U and J. We 
find that the main quantity that determines the gap is 
the inter-orbital Coulomb interaction U' = U — 2 J, and 
we get A « 1.2 eV for U = 3.3 eV, J = 0.5 eV. In Fig. [3] 
we plot the spectral function for different sets of inter- 
action parameters, and find A r* 1.4 eV for U = 4eV, 
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FIG. 4: (Color online) Imaginary part of the self-energies on a 
reference system consisting of two coupled 4-site chains. Left 
panel: Real space £ij(a;). Solid line: Local/onsite self-energy 
£n(u>). Dashed line: Inter-chain self-energy £12 (w). Dash- 
dotted line: Intra-chain self-energy £13 (tj). Right panel: Self- 
energy for selected momenta of the cluster. Momenta are in- 
dicated in the plot, and a vertical shift between momenta has 
been introduced for improved presentation. The self-energy 
shows causality (negative definite), almost no dependence in 
a-direction, and strong dependence in o-direction. Lorentzian 
broadening of 77 = 0.02 eV was used in both plots. 



J = 0.7 cV, and A w 1.9 cV for U = 4cV, J = 0.5 eV. 
All these values for the gap are a bit smaller than the ex- 
perimental charge gap of about 2 eV^ but nevertheless 
in reasonable agreement. 



B. 



Spectral weights in the single-band Hubbard 
model 



We have shown in Sec. IIII Al that the occupied states 
of the t2 g manifold are well reproduced by a single-band 
Hubbard model. In this section we want to investi- 
gate the spectral function of Eq. Q in more detail. As 
mentioned above, we focus on the effect of the addi- 
tional two-dimensional hopping processes on the quasi- 
one-dimensional behavior of TiOCl. 

First, we want to determine the strength of the correla- 
tions along the qualitatively different bonds of the lattice, 
Fig. m This can be done best by inspecting the self- 
energy Tiij(u>), which is, in VCA, a quantity defined on 
the reference system and, thus, can be readily obtained. 
The results for three selected matrix elements arc shown 
in Fig.[31 calculated on an 2 x 4 cluster. This cluster con- 
sists of two coupled 4-site chains in 6-direction. In other 
words, a cluster similar to the one depicted in Fig. [T]i) 
but with doubled extension in 6-direction. It is obvious 
that local (En(w)) and intra-chain correlations (£13(0;)) 
are strong, but the correlations between adjacent chains 
(£12 (w)) are orders of magnitude weaker. In addition, 
we show in the right panel of Fig. [3] the Fourier transfor- 
mation of the self-energy, for momenta accessible at this 
small cluster. Again, the influence of momenta perpen- 
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FIG. 5: Spectral function of the single-band model with a 
12-site chain as reference system. Top: No coupling between 
chains. Bottom: Inclusion of interchain coupling parameters 
obtained from LDA-downfolding. Dark areas mark large spec- 
tral weight. Lorentzian broadening of 77 = 0.02 eV was used. 
The horizontal striped structures occur due to the breaking 
of translational invariance in the cluster approximation. 



dicular to the chains is hardly visible in the self-energy, 
whereas it shows significant momentum dependence in 
the chain direction. This leads to the conclusion that the 
spectra of TiOCl should be governed by ID correlations, 
modified by single-particle effects due to the coupling of 
the chains. 

Motivated by this result, we use from now on a 1 x 12 
cluster as reference system. Since the VCA approx- 
imates the interacting Green's function as G(w) _1 = 
Gpj - E(w) with G , t the non-interacting Green's func- 
tion of the model Hamiltonian and the cluster self- 
energy, it is easy to see that the inter-cluster coupling 
is treated in a single-particle (i. e., non-interacting) man- 
ner, since it enters just via Go,t- On the other hand, 
this procedure gives the best possible description of the 
correlation effects along the chains in 6-dircction. 

In Fig. [5] we show a density plot of the spectral func- 
tion of the Hamiltonian Eq. (J3j> for U = 3.3 eV. In the 
upper part we included only the intra-chain hopping 
ti = — 0.21eV in the calculation, leading to flat bands 
in a-direction, i. e., from (vr, 0) to (0, 0), since in this case 
the chains are decoupled. By including additional inter- 
chain parameters ti = 0.03 eV and t% = 0.04 eV as given 
in Ref . [llil. we notice a slight dispersion in a-direction. In 
6-direction, however, the band positions remain almost 
unchanged; we find only changes in the spectral weights. 
Since this cannot be seen clearly in the density plots, 
we show in Fig. [5] the evolution of the spectral function 
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FIG. 6: Spectral function A(k = 0, to) at the P point. Top left: 
without inter-chain coupling. Top right: Including nearest- 
neighbor inter-chain coupling t2- Bottom right: Including 
next-nearest- neighbor inter-chain coupling ts,. For compar- 
ison, bottom left: next-nearest-neighbor hopping along the 
chain, but no inter-chain hopping. Lorentzian broadening of 
n — 0.02 eV was used. 



at the r point k = (0, 0) when longer-ranged hopping 
processes are included. 

The upper left panel (a) is the spectral function for 
decoupled chains with the spin-charge separation clearly 
visible. At the T point the holon band is located around 
— 1.72eV, and the spinon band around — 1.39cV. Includ- 
ing the nearest neighbor inter-chain hopping t2 leads to 
a significant redistribution of spectral weight from the 
holon to the spinon band, i. e., from higher to lower bind- 
ing energies, see the upper right panel (b). This effect is 
even enhanced when the next-nearest inter-chain hopping 
£3 is included, as shown in the lower right panel (d). Here 
the spectral weight of the low binding energy ('spinon') 
excitation is comparable to the weight of the high binding 
energy ('holon') excitation for decoupled chains in panel 
(a), and vice verse. At this point we would like to men- 
tion, that in a strict sense the terminology 'spinon' and 
'holon' is not applicable any more, since these are prop- 
erties of purely onc-dimcnsional systems. Anyway, since 
the spectra resemble to some extent ID systems, we still 
use these terms to distinguish the different excitations. 

From Fig. [5] it is clear that the inclusion of inter-chain 
processes enhances the asymmetry of the k-resolved spec- 
tra. The low-lying excitation near the T point is strongly 
enhanced, whereas there is no spectral weight transfer to 
lower binding energies visible around (0,7r). Note that 
we define the spectrum to be symmetric if the main ex- 
citations at k-vector (0, 0) and (0, 7r) are located at the 
same binding energy. 

One may ask if it is possible to get a similar spec- 
tral weight distribution by using only the purely one- 
dimensional Hubbard model, but including longer-ranged 
intra-chain hopping processes as given by the ab initio 




(0.0) (0,lt) (0,0) (0,n) 



FIG. 7: (Color online) Comparison of the theoretical spectral 
function (left) with the experimental ARPES spectra (right). 
Only the ^-direction is shown. Lorentzian broadening of n = 
0.02 eV was used. The vertical dashed line marks momentum 
(0, 7r/2), and the band widths wt and w' b are indicated. Dark 
areas mark large spectral weight, the color scale is normalized 
in each plot separately. 



calculations. In fact, the next-nearest-neighbor hopping 
term along the chain, £4, is of similar size of the inter- 
chain hoppingSf^ The result for the spectral function at 
the r point in this pure ID case is shown in the lower 
left panel (c) of Fig. [5] From this result it is obvious 
that one can not get an excitation at binding energies 
of roughly — 1.4eV, as seen in experiments. On the con- 
trary, the spectral weight of the excitation at the higher 
binding energy of about — 1.8eV is even enhanced in the 
one-dimensional treatment when longer-ranged hopping 
processes are included. 

Our results support the description of TiOCl as a lay- 
ered two-dimensional compound with strong anisotropy, 
also on the level of correlations. There is finite disper- 
sion also perpendicular to the chains, but a backfolding 
of the bands can only be seen along the chains, where 
correlations are dominant. 

Let us now comment on the relation of our results 
to ARPES data. Experiments showi^ that the disper- 
sion in TiOCl shows a strong asymmetric behavior along 
the crystallographic b direction, see also the right plot 
of Fig. [7J The binding energy of the lowest-lying band 
around k = (0,0) is about — 1.5eV, whereas around 
k = (0,7r) it is about — 2.0eV. First attempts to de- 
scribe the dispersions within ab initio calculations were 
not successful. Standard LDA calculations do not pro- 
duce the backfolding of the bands induced by short-range 
spin fluctuations, and spin-polarized LSDA+U calcula- 
tions cannot account for the asymmetry of the spectra. 
Also the spectra of the one-dimensional Hubbard model 
calculated within the dynamical density-matrix renor- 
malization group (DDMRG) do not reproduce the exper- 
imental spectral weight distribution^ However, as can 
be seen in Fig. our new results show that the Hubbard 
model can indeed give a good description of the asymme- 
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try. since the inter-chain processes give rise to a spectral 
weight transfer from the holon to the spinon band around 
the r point, and therefore make the band structure more 
asymmetric. 

In order to compare our results more quantitatively, we 
extract the following numbers related to the band widths 
of the spectra. The first one, Wf,, is the difference of the 
binding energies at (0,7r/2) and (0,7r) and is a measure 
for the overall band width. The second one, w' b , is de- 
fined to be the difference in binding energies between 
(0, 7r/2) and (0,0). The larger the difference between 
these two quantities, the larger is the asymmetry of the 
spectra. From experiments^ 5 - we extract Wb ~ 0.47 eV 
and w' b f» 0.17cV, and for the calculated spectra we get 
Wb ~ 0.50 cV and w' b ~ 0.09 eV. Comparing theory and 
experiment, we see that the overall band width Wb is 
well reproduced by the calculation, but the asymmetry 
is even a bit more pronounced in the theoretical spec- 
tra, resulting in a smaller value of w' b . This result may 
be improved by including more longer-ranged hopping 
processes. Although they decrease rapidly with distance, 
they can change the band widths within a few percent. 

By inspecting Fig. [7] it is obvious that the width of 
the spectra is much larger in the ARPES data than in 
the calculated spectra. The most important reason for 
that is that our calculations are done at T = 0, using 
exact diagonalization techniques. Therefore there are no 
life-time effects due to finite temperatures included in 
our calculations. Moreover, additional coupling to lattice 
degrees of freedom could also lead to a smaller life time, 
and hence broader excitations. 

In summary, our results imply that it is important to 
include the inter-chain couplings at least on a single- 
particle non-interacting level into the effective model, 
in order to improve the description of the experimental 
spectra, although the strong correlations are constricted 
mainly to the chains. 

With our work we could fill the gap left by previ- 
ous theoretical studies regarding the momentum-resolved 
single-particle spectral function of TiOCl. There are, 
though, still some open questions. For example, we do 
see clear signatures of spin-charge separation in our calcu- 
lated spectra, which have not been found experimentally. 
Also the so-called shadow band, dispersing at around 
—2 eV has not been seen in the ARPES spectra. A reason 
for this can be a very small relative spectral weight of the 
high-energy band that cannot be resolved in experiment. 
In our calculation we also did not include the lattice de- 
grees of freedom, which are supposed to be very impor- 
tant in TiOClii, driving the spin-Pcicrls phase transition. 
These phonons can also lead to a smearing of the peak 
structure of A(k, u>). 

Finally, we want to comment also on the differences be- 
tween the ARPES spectra of TiOCl and TiOBr. Exper- 
iments, supplemented with band structure calculations, 
have show n 15 ' 45 that in the latter compound the intra- 
chain couplings are weaker and the inter-chain couplings 
stronger compared to TiOCl^ e.g., t\ decreases from 



— 0.2IcV to — O.I7cV, whereas t% increases from 0.04eV 
to 0.06 eV. By inspecting the self-energies in a similar 
manner as we did in Fig. [H we found that also for these 
parameter values the correlations are predominantly one- 
dimensional. There are only changes in the overall band- 
widths, but no qualitative changes. For instance, the 
band width in a-direction is enhanced, but there are no 
signatures of strong inter-chain correlations resulting in 
a backfolding of the bands. 

At this point we want to remark that, in particular 
for TiOBr, one should be very careful with the use of 
the spinon/holon terminology. Although there are no 
qualitative changes due to the enhanced couplings, they 
do change the band widths. Hence, quantitative analysis 
have to be done including these inter-chain couplings. 

IV. CONCLUSIONS 

By combining ab initio calculations (LDA) and the 
variational cluster approximation, we could study for the 
first time the momentum-resolved spectral function in- 
cluding a realistic band structure and strong-correlation 
effects. In agreement with previous theoretical studies 
and experimental results, our calculations showed an al- 
most complete polarization in the orbital sector, with 
99% of the electrons occupying the Ti d xy orbital. 

Since the orbital degree of freedom is quenched, we 
could use an effective single-band model for the investi- 
gation of the spectral properties. The most striking result 
of our study is that the inclusion of inter-chain hopping 
processes leads to a significant spectral weight redistri- 
bution around the T point from higher to lower binding 
energies. This effect, which makes the spectrum more 
asymmetric with respect to the points (0,0) and (0, 7r), 
cannot be reproduced using only the hopping processes 
along the chains. This result suggests that the frustrated 
inter-chain coupling^ is one of the main reasons for the 
strong asymmetry that has been found in experimental 
ARPES measurements. Moreover the calculated spec- 
tral features may be extended to a more general class of 
correlated coupled-chains systems. 

An open question in TiOCl is still the role of phonons. 
Because of the vicinity of the system to a spin-Pcierls 
instability, the phonons are supposed to be important in 
the system. Including these degrees of freedom, although 
theoretically very demanding, could further improve the 
results with respect to lineshapes and lifetimes of the 
excitations. 
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